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Three interrelated questions concerning Kerr spacetime late-time scalar-field tails are considered 
numerically, specifically the evolutions of generic and non-generic initial data sets, the excitation of 
"up" modes, and the resolution of an apparent paradox related to the superposition principle. We 
propose to generalize the Barack-Ori formula for the decay rate of any tail multipole given a generic 
initial data set, to the contribution of any initial multipole mode. Our proposal leads to a much 
simpler expression for the late-time power law index. Specifically, we propose that the late-time 
\ decay rate of the Yim spherical harmonic multipole moment because of an initial Yfm multipole is 

fSJ ■ independent of the azimuthal number m, and is given by where n = t' + I -\- 1 for I < t! and 

n = ^' + ^ + 3 for £ > £' . We also show explicitly that the angular symmetry group of a multipole 
O ' does not determine its late-time decay rate. 
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I. INTRODUCTION AND SUMMARY 



Perturbations of black holes decay first with complex frequencies, known as the quasinormal modes of of the black 
5—1 hole, taken over at late times by the tails of the perturbation field — decaying as power-laws of time. The decay 
rate of the late-time tails in the spacetime of spinning black holes has been the subject of much debate and some 
confusion. Much of the debate in the literature was focused on the late-time decay rate of an initially pure azimuthal 
CNJ ■ hexadecapole {£' — 4, m — 0) scalar field perturbation. All authors agree that at late times the field is dominated 
^ ] by its monopole moment, and that the decay rate would be according to an inverse power of time, but different 
' , ^ i: claims were made as to the value of the power-law exponent. Specifically, there were claims for decay rate along fixed 
^ ' Boyer-Lindquist r according to Q, [3,I3|, and even t~^-^ 

In the last couple of years the behavior of the power-law tails has been clarified in a number of papers 

143 (and 

the reader is referred to the detail therein), yet there are some remaining interesting detail waiting to be unveiled. 
In particular, it was shown in f?\ that the disagreement between the and behaviors can be attributed to the 
use of different initial data sets (see also but not to the use of different slicing conditions as was suggested in 
0. Most importantly, Boyer-Lindquist slicing is not unique in terms of the resulting tails as had been argued, but 
rather belongs to an equivalency class of slicing conditions (including also ingoing Kerr slicing), all having similar tail 
structure properties. 

For the evolution of scalar fields on a fixed background of a Kerr black hole, the value of the power-law index, 
namely the exponent n in the late-time decay rate was given by Hod Q as n = £' -|- m -f 1 for even £' — m > 2, 
5^ I n=£' + m + 2 for odd £' — m > 2, and n = 2£' + 3 for £' — m < 2, where £' is the multipole moment of the initial 
perturbation and m is its azimuthal number. This expression for n is rather complicated and notably depends on the 
azimuthal number m. We argue that the added tacit requirement that it is the decay rate of the slowest decaying 
mode that is of usual interest masks the inherent simplicity of n in the Hod formula. Specifically, we propose that 
the decay rate of the £ mode because of an initial £' mode is independent of to, and is given by rt = -I- ^ -I- 1 for 
£ < £! and n = £' + £ + for £ > £' . Notice that this proposal is consistent with Hod's formula and with when 
they are applicable, and generalizes them also to all other cases. It is therefore a proposal for the general formula for 
the decay rate of any multipole moment of the tail because of any multipole moment of the initial data. It is striking 
that this more general formula is simpler than the formula for the decay rate of the least damped mode, that was 
the question that 0] was addressing. We believe that this simple, general formula will provide insight into a more 
complete understanding of the tail phenomenon in spinning black hole spacetimes. 

Our general formula allows up to address some questions left open also in the analysis of Barack and Ori Q. Barack 
and Ori studied analytically the late-time tails for generic families of initial data (only assuming compactly supported 
outgoing initial pulses) . Instead of carrying the analysis in the frequency domain — which turns out to be complicated 
on a Kerr background because the frequency dependence of the spheroidal harmonics implies that separation of the 
two angular variables depends on the frequency — Barack and Ori analyzed the evolution of perturbations in the time 
domain. Barack and Ori found the decay rate of the £ multipole given generic initial data, which they considered to be 
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data in which all multipoles are present. In fact, because the smallest multipole contribution to the initial data turns 
out as we show below to dominate the decay rate of the £ multipole at late times, Barack and Ori found the decay 
rate of an even (odd) mode i excites by the scalar field monopole (dipole) mode of the initial data. The Barack-Ori 
formula for the decay rate of the t multipole has not been confirmed by numerical simulations, which we do below for 
the first time. Our proposed formula also generalizes it to any initial multipole, not just the monopole (or dipole) as 
in 0. 

Consider next the question of the decay rate of a certain multipole moment I. In the Schwarzschild spacetime 
the decay rate depends only the the value of specifically it is where n = 2£ + 3. Most notably, the decay 
rate is independent of the history of the mode. This is no longer the case in the Kerr spacetime: the decay rate 
of a multipole I does depend on how that mode came into existence, specifically, it depends on both i^l' . In the 
Schwarzschild case the dependence is only on i because of the degeneracy I' — I and the absence of mode couplings 
("spherical harmonics are eigenfunctions of the wave operator") . Therefore, in the Schwarzschild case one may say 
that "a multipole is a multipole," and consequently its decay rate is intimately linked to its angular distribution and 
the associated symmetry group. Specifically, in Schwarzschild there is a one-to-one relation of the symmetry group 
of the mode in question and the mode's late-time decay rate, such that the symmetry determines the time evolution 
of the mode. Previous studies of tail behavior in Kerr have suggested that it is not correct to say in Kerr that "a 
multipole is a multipole," although the authors are not aware of explicit statements in this regard. We bring here 
new evidence that that is indeed the case, and the tail decay rate has nothing to do with the symmetry group of the 
mode in question. This insight may be helpful in the understanding of cases where more than one channel contributes 
to a mode of a certain i value, although each channel contributes so that the decay rate is different, in what might 
be construed as an apparent violation of the superposition principle. The realization that in Kerr one may no longer 
correctly argue that "a multipole is a multipole" completely resolves this apparent violation. 

The organization of this Paper is as follows: In Section|lI]we describe the numerical code and the type of initial data 
sets that we use. Specifically, we describe the code development that was necessary in order to support our proposal 
for the tail decay-rate formula. In this Paper we discuss three related questions: First, in Section HII] we address the 
difference in the time evolution of late-time tails between generic and non-generic initial data sets. Generic initial 
data sets were first considered by Barack and Ori l8|, but their results have not been confronted with numerical 
simulations. In addition to verifying the results of Q, we also discuss the meaning of genericity of initial data sets. 
As we demonstrate, not any mixture of all multipole modes is in fact generic, and one needs to exclude mixtures 
that are the outcome of the evolution of non-generic initial data, although such mixtures exhibit similar properties to 
mixtures that represent hona fide generic data sets. Then, in Section ITVl we consider an aspect of late-time tails that 
has not been studied before, specifically the decay rate of excited higher multipole modes, which we call "up" -excited 
modes, namely the case £ > i' . The decay rate of "up" excited modes generalized the Barack-Ori result in Q to 
some interesting cases: the case of non-generic initial data in which (only) the lowest dynamically allowed mode is 
not present, and the case in which only one multipole mode is present in the initial data. 



II. NUMERICAL CODE 



A. Initial data sets 



Linearized perturbations of Kerr black holes are described by Teukolsky's master equation, given in Boyer-Lindquist 
coordinates for a scalar field (s = 0) ^, known as the Teukolsky function, by 
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where M, a are the black hole's mass and spin angular momentum per unit mass, respectively, and the horizon function 



A = 



2Mr + a'^. 



To solve Eq. ([T]) numerically, we define, following Q, 5(r, 6) :— {r^ + a^)/I] where = (r^ + a^)^ — A sin^ I 
We next define the 'momentum' 11 of the field <j) defined hy tp — e""''' (p, according to 
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where r, is the regular Kerr spacetime 'tortoise' coordinate defined by dr^j dr = (r^ + a^)/ A. 

Our code is a 2+lD first-order code for the time evolution of modes to of (/), 11 given two independent functions of 
two variables at i = as initial data, based on Q. The convergence order is taken to be second-order in the radial 
and temporal directions, and sixth order in the angular direction (although some of the higher-^ results were obtained 
with tenth-order angular operators). The standard grid resolution we use is 64, 000 grid points in the radial direction, 
and 64 angular grid points (which we denote as QAK x 64), unless stated otherwise, taking the temporal step to be 
the largest step consistent with the Courant condition. 

In practice, we choose as initial data (j>{r,9)\t=o = and Il{r,9)\t=o — /(r) ( cos 0), where the radial function 
/(r) is chosen to be a gaussian /(r) = (I/v^tht^) exp[— (r, — o)^/(2o'^)] with o = 25M and cr — 6M. The outer 
and inner boundaries are placed at = ±800Af, which allows us to integrate to t = 1, 500 Af at r* o without 

seeing boundary refiection effects. We choose the Kerr parameter a = 0.995Af. We use quadrupole or octal precision 
floating-point arithmetic throughout according to need. 

We emphasize that our choice of initial data — based on the field and its conjugate momentum — is in accord with 
other common choices. Specifically, when we choose a single value of i, our choice of initial data is a pure multipole 
moment, equivalent to choices made based on the field and its time derivative. There is no ambiguity in the definition 
of a pure multipole [7| . 



B. Code Improvements 

Since the goal of this work is to investigate the late-time decay rate of higher multipole modes of Kerr black holes 
a number of improvements to our time-domain Teukolsky equation code were necessary. This is because generating 
accurate numerical data in the context of higher multipole modes involves a number of challenges. Firstly, these 
simulations need to be rather long — this is because typically the observed field exhibits an exponentially decaying 
oscillatory behavior ( "quasi- normal ringing" ) in the initial part of the evolution and only much later this transitions 
over to a clean power-law decay. One needs to wait for these oscillations to dissipate away. For higher multipoles 
this phase lasts much longer, thus requiring significantly longer evolutions. Secondly, because each multipole has a 
different decay rate (which increases with an increase in £) at late times one ends up with numerical data in which 
different multipoles have widely different amplitudes (often 30 - 40 orders of magnitude apart!). For this reason, not 
only does the numerical solution scheme have to be high-order (to reduce the discretization errors to the required 
levels) but it also requires high-precision floating-point numerical computation (due to the large range of amplitudes 
involved) . 

Below we list the major improvements made since the publication of [7] to our Teukolsky equation evolution code: 

1. Higher-Order Algorithm: As mentioned above, we require a higher-order numerical evolution scheme to keep 
the discretization errors to sufficiently low levels, otherwise these can overwhelm the numerical data associated 
to the weakly excited higher multipole modes of interest. It turns out that it is sufficient that only the angular 
differentiation (i.e. ^-derivatives of the field) be implemented using a higher-order numerical stencil. The 
temporal and the radial direction related operations can simply stay 2nd-order and such a mixed approach yields 
sufficiently good results For this reason, in this work we choose the finite-difference angular differentiation 
operator to be lOth-order accurate and leave the rest of the numerical scheme as a standard 2nd-order Lax- 
Wendroff algorithm. 

2. Higher- Order Projection Integrals: We also require a matching lOth-order accuracy in mode projection integrals 
to prevent numerical integration errors from overwhelming the higher multipole numerical data. 

3. Higher-Order Numerical Precision: As pointed out above, we also require high-numerical precision in our 
computations. In particular, double, quadruple and octal precision may be required depending upon the details 
of Kerr tails simulation being attempted. This is simply to prevent roundoff error from overwhelming the 
weakly excited modes of interest - especially important for "up" excited higher multiple modes! We made use of 
quadruple precision (in some cases, octal precision) floating-point accuracy in this work. Now, very few current 
computer processors support full quadruple-precision (128-bit) datatype and operations. And to the best of our 
knowledge, no compute hardware natively supports octal-precision (256-bit) arithmetic. Thus, we implemented 
such high-precision computations in software using publicly available software libraries. This made our work 
particularly challenging computationally, because with each such increase in numerical precision, our code's 
overall performance drops by an order-of-magnitude. 

4. Code Parallelism and HPC Hardware: In order to perform the numerical simulations in a reasonable time- 
frame, we made use of message-passing (MPI) based parallelism in our code and also used sophisticated HPC 
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hardware. In particular, domain-decomposition based parallelism was implemented in the code and IBM Cell 
BE and Power? systems with specialized optimizations fld\ were used to perform these long evolutions in a 
reasonable time-frame. 

5. Other Enhancements: We also enhanced the code further to be able to evolve non-axisymmetric mode compu- 
tations (both odd and even m-modes). The even m-mode case was straightforward to implement, since these 
simply require some additional terms in the evolution equation. For the odd m-mode case we used the follow- 
ing procedure to simplify our implementation: Instead of evolving the scalar field ip itself, we transformed the 
Teukolsky equation to evolve the field tp/ sin 9. This simple transformation enabled us to reuse several of our 
even m-mode code algorithms (such as the higher-order angular boundary condition imposition) and also avoid 
the introduction of non-polynomial functions of the code's angular grid variable, i.e., cosO. 

III. NUMERICAL EVOLUTION OF GENERIC INITIAL DATA 

We first show for the first time numerical tests of the Barack-Ori formula 8] . Then, in Section IIVI we show how 
the Barack-Ori formula is a particular case for "up" excitation of modes. Barack and Ori noted that in the generic 
case of black hole perturbations, all dynamically allowed modes are present. Specifically, in the case of scalar field 
perturbations, all modes, starting with the monopole {£ — 0) evolve. The decay rate of a multipole i with azimuthal 
number m is given by t""' where 

ne = e+ |m| + 3 + g , (2) 

where q — 0,1 ii £ + m is even or odd, respectively. Notably, the Barack-Ori formula depends explicitly on the 
azimuthal number m, and like the Hod formula includes a certain complexity that has to do with whether £ + m is 
even or odd. 

The meaning of "generic" initial data is that the relative amplitudes of the various multipoles and their radial 
profiles are not fine tuned. Consider, e.g., the following initial data set, 0(r, 0)|t=o = and n(r, 0)t=o — A sin^ 9 f{r), 
where /(r) is any localized radial function, say a gaussian. These initial data include non- vanishing monopole and 
quadrapole moments. We evolved these data numerically, and found that as expected the late-time tail includes all 
even multipoles, and is dominated by the monopole field that decays according to t~^ at late times. The quadrapole 
projection of the late-time field decays according to t~^, in accordance with the Barack-Ori formula. Other choices of 
initial data sets and multipoles are also found to be in agreement with the Barack-Ori formula. As far as the present 
authors are aware, this is the first direct numerical confirmation of the Barack-Ori formula. 

Let us now consider the following scenario: one starts with initial data as above, and evolves the field to a certain 
value of time, say to t = T. The field at t = T is a complicated mixture of modes, each with its own radial profile. 
Now consider a different initial value problem, in which at i = only a pure multipole, say the quadrapole moment, 
is non-vanishing, and evolve these initial data to t — T. These pure mode initial data also evolve to a complicated 
mixture of multipole moments at i = T, and at first look one does not notice much qualitative difference between 
the fields at t = T. However, treating the fields at t — T as new initial value problems, further time evolution leads 
to different decay rates at late times. Specifically, the generic data lead to quadrapole projection that drops off like 
t~^, and the pure data lead to quadrapole projection that drops off like t~^. Considering only the values at t = T as 
the initial data sets, what is the fundamental difference between the two sets, that leads to very different late-time 
evolutions? Even though either set appears at i = T to be a set of amplitudes with weighted ratios for the various 
multipole moments (as functions of the radial coordinate), the set describing the pure mode evolution is made at 
t = r of a very carefully chosen set of (r dependent) amplitudes. In fact, it is the outcome of fine-tuned evolution 
of a pure mode of the original initial value problem. Examination of the data set at t = T does not reveal anything 
qualitative different about it: it is the specific amplitude ratios and radial profiles that make it non-generic. 

One may test these ideas by injecting at t ~T another field, say that of a pure mode. Start at t = Q with a pure 
quadrapolar field. Kt t = T the projection of the field is a mixture of all even multipoles. Then, at t — T, inject either 
a pure monopole or a pure quadrapole with some radial profile. The multipolar content of the field after the injection 
is changed in that the relative multipolar amplitudes changed (as functions of the radius). When the injected field is 
quadrapolar, the late-time decay rate of the quadrapole moment is as should be expected from the fact that we 
have a superposition of two linear pure mode evolutions (which a certain time delay between them) (Fig. [1]). When 
the injected field is a monopole, the late-time decay rate of the quadrapole moment is t^^ (Fig. [2]). The addition of 
a monopole field (without fine tuning its amplitude and radial profile) makes the multipolar content of the total field 
generic. One may pose the following question: When the injected field is a quadrapole, its time evolution will excite 
a monopole. Why in the presence of this excited monopole the quadrapole projection of the field decays like t~^, 
whereas in the presence of an injected monopole the decay rate is the slower t~^? There is a fundamental difference 
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FIG. 1: The quadrapole projection of the full field as a function of time, for initial data of a pure quadrapole field, then injected 
at T = 250M with a pure quadrapole field. The late-time field behavior is computationally found to be given hy ip ^ t~^'^^. 
The power law index, and all the indices in this Paper, were calculated by extrapolating the local power-law index to infinite 
time. 




FIG. 2: The quadrapole projection of the full field as a function of time, for initial data of a pure quadrapole field, then injected 
at T = 250M with a pure monopole field. The late-time field behavior is computationally found to be given by ~ 

between the two eases: the exeited monopole is not an arbitrary field: it is earefully chosen by the dynamies of the 
problem, and may be viewed as a fine tuned field; the injected monopole is arbitrary, and is not fine tuned to lead to 
a different late-time decay rate. 

IV. EXCITATION OF "UP" MODES 

Most interest has been devoted to finding the decay rate of the slowest decaying mode, for obvious reasons. Specif- 
ically, the slowest decaying mode determines the late-time behavior of the full field. As lower multipole modes 
typically have a slower decay rate, most interest has naturally focused on the excitation of "down" modes, specifically 
the excitation of lower multipole modes starting with higher multipole modes. 

Here, we are interested in the converse, i.e., in the excitation and decay rate of "up" modes, i.e., the excitation and 
decay rate of a higher multipole mode than the one initially excited. Such modes are typically sub-dominant, and do 
not dominate the late-time behavior of the full field. Our interest in the "up" modes is therefore only tangentially 
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FIG. 3: The I = 0, 2, 4 multipole projections, for initial data of a pure quadrupole field as functions of time. The late-time 
field behavior is computationally found to be given by 2'^o ~ ^-^-^^^ ^ 2ip2 ~ £~^"°^, and 24'i ~ t~^'°^, respectively. 

related to the question of the decay rate of the full field. Instead, it is the gaining of more insight towards a more 
complete understanding of the dynamics and mode coupling properties thereof we are interested in here. Denoting 
the initially excited pure multipole by i' , and the multipole moment of the field whose excitation and decay are is of 
interest by £ (so that for "up" -modes £ > £' and for "down" -modes £<£'), we propose that the late-time decay rate 
for excited £, m modes along a r = const worldline, is given by t~" where 

_ ( £' + £+l foT £<£' , . 

+ ^ + 3 for £>£' ■ 

Notably, this formula is independent of m. Specifically, for all (allowed) values of m the decay rate of the £ excited 
multipole because of an initial £' is the same. This formula is simpler than either the Hod or the Barack-Ori 
formulas, but is fully consistent with them in their domain of applicability. Specifically, when £' — 0,1, Eq. (|3]) 
coincides with the Barack-Ori formula ([2]). In the Schwarzschild case there is no mode coupling, £ = £' , and therefore 
n = £'+£ + i = 2£ + i. The numerical evidence described below supports Eq. ([3]) in all cases we have checked. 

Azimuthal modes 

To test our proposal — and in particular the 771-independence of the decay rate, which may be construed as a 
statement on the generally m-dependent Green's function — we focus now on pure multipole initial data sets. 

Letting the initial £' = 2, we project the multipoles ^ = 0,2, 4, and the local power indices are shown in Fig. [3l Of 
particular interest here is the excitation of the azimuthal, hexadecapole mode, £ = A. In this case, £' = 2,£ = 4, so 
that our proposal is that n = 2 + 4 + 3 = 9. In our numerical simulations, we have found the value of — 9.01, 
with uncertainly in the last figure. 

Starting with initial data for a pure £' = {) monopole mode, we project the £ — Q,2 and 4 modes. The local power 
indices o«o and 0^12 are shown in Fig. H) The hexadecapole mode £ — 4, 0^4, is shown in Fig. [S] Figure [5] suggests 
that our grid density is not sufficiently high to obtain accurately the very late time hexadecapole projection to the 
required accuracy. Indeed, as Fig. [S] shows, the tail of the field "curves up" at very late times. However, increasing 
the grid density appears to straighten up the field, from which we infer that this curving is a numerical artifact. 
We attribute the problem to our calculation of a mode that is a second-order excited "up" mode, i.e., the £ = A 
excited mode from an initial £' = mode. The practical problems in determining the decay rate accurately are 
accentuated by Fig. |6l which displays the local power indices for three grid densities. Figure [6] suggests that indeed 
a higher resolution simulation would be successful at determining the late-time decay rate with sufficient accuracy. 
Our results do suggest that qu^ — 7 (notice, this is a numerical result with one significant figure, not an exact value!), 
as indicated by the reference curve in Fig. [5l but any determination of the decay rate at higher accuracy than that 
would require higher resolution simulations. 

We have also tested our proposal with odd modes. Specifically, Let £' = 3 and study the late time tail for the modes 
£ = 1,3, 5. Figure [7] shows the projections ^ = 1, 3, 5 as functions of time, and figure [S] displays the local power indices 
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FIG. 4: The local power indices o^o (upper panel) and 0^2 (lower panel) as functions of the time t. 




FIG. 5: The £ = 4 multipole projections, for initial data of a pure monopole field as functions of time. Shown are the fields 
with three different grid densities, 64K x 64 (dash-dotted curve), 80K x 80 (dashed curve), and WOK x 100 (solid curve). The 
inserts show the same for segments of the full data, in addition to a reference curve ~ t~'^ . 



TABLE I: Scattering matrix for azimuthal modes. Late-time power law indices i'Ue for the I'-projections of the fields for pure 
azimuthal (m = 0) mode initial data sets, £' — 0, 1, 2, 3, 4. This table includes both "down" and "up" excitations. Asterisks 
relate to kinematically disallowed states. All figures are significant. 
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FIG. 6: The local power index on4 as a function of the time t for three grid resolutions, 64:K x 64 (dash-dotted curve), 80K x 80 
(dashed curve), and WOK x 100 (solid curve). 



5 




FIG. 7: The ^ = 1,3, 5 multipole projections, for initial data of a pure octupole field as functions of time. The late-time field 
behavior is computationally found to be given by ipi ~ f"^""^, ip^ ~ ^-9.009^ ^ ^-ii.oos^ respectively. 



for these modes as functions of time. 

Starting with £' ~ 1, we project the dipole and octupole modes and their respective local power indices in Fig. |9l 
We encounter the same problem finding the decay rate for £ — 5 — i.e., find iip5 and in^ — as we did when we 
we calculated oipi and o«4 (Figs. fTOl and [TT|) . Again, the second-order "up" -excited mode requires high resolution 
simulations. 

Our proposal ([3]) naturally extends the Barack-Ori formula ([2]) to cases where the initial data do not include the 
monopole {£' — 0) mode. Recall that the Barack-Ori formula considers the case of generic initial data, which means 
a mixture of all multipole modes with arbitrary relative amplitudes and radial profiles. In particular, the monopole 
mode is generically present, so that the slowest decaying contribution of any "up" -excited £ mode is generated by 
the initial monopole. Our proposal degenerates to the Barack-Ori formula when £' = 0. When the monopole is not 
present in the generic initial data set, our proposal will produce the decay rate of a mixture of modes when £' is taken 
to be the smallest multipole present in the generic data set, and as a particular case will also produce the decay rate 
of the £ multipole when the initial data is a pure £'{< £) mode, the non-generic case, in the language of Barack and 
Ori. In summary, the Barack-Ori formula Eq. ([2]) is simply the particular case on£ of Eq. ([3l). 

Table I summarizes our results for the azimuthal, m = 0, case. It shows the power law indices for the late 
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FIG. 8: The local power indices of the £ — 1,3,5 multipole projections as functions of time, for initial data of a pure octupole 
field (same data as in Fig. [7} • 
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FIG. 9: The dipole {£ = 1) and octupole {£ = 3) projections, i-i/^s and iips, correspondingly (upper panel), and their local power 
indices, ini and ina, respectively (left and right lower panels), as functions of time, for initial data of a pure dipole {£' = 1) 
field. 

time tails for both "up"- and "down" -modes. It is in some loose sense a "partial scattering matrix," that includes 
(partial) information on the excitation amplitude for each of the excitation channels corresponding with a preparation 
state of the initial data. The scattering matrix is only partial, because it includes no information on the pre-factor 
"00 = V'o(''Oj profile) of the late-time tail ip(t) ^ ipo i^", where tq is the evaluation point, and where ipo also possibly 
depends on the radial profile of the initial data. The (partial) scattering matrix includes only the information on n. 
The scattering matrix has a clear shape of the sum of two triangular matrices: letting rows be associated with £' and 
column with £, the upper triangular matrix elements satisfy n — £' + £ + 3 and the strictly lower triangular matrix 
elements satisfy n = £' + £+ 1. 

Non— azimuthal modes 

Our proposal ([3]) is not limited to azimuthal modes, and predicts that the late-time decay rate of the £, m modes 
that is excited by the initial £', m mode is independent of m, and depends only on t, I' . We have tested this hypothesis 
in numerous cases, that are summarized in Tables I, II, and III. 
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FIG. 10: The £ — 5 multipole projections, for initial data of a pure monopole field as functions of time. Shown are the fields 
with three different grid densities, GAK x 64 (dash-dotted curve), 80K x 80 (dashed curve), and lOOA" x 100 (solid curve). The 
inserts show the same for segments of the full data, in addition to a reference curve ~ . 
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FIG. 11: The local power index ins as a function of the time t for three grid resolutions, GiK x 64 (dash-dotted curve), 
80K X 80 (dashed curve), and lOOA' x 100 (solid curve). 



In Fig. [12] we show the fields for m = I and m = 2 modes that are created by an initial octupole and excites all 
odd £ modes. Specifically, we show the £ — 3,5 projections for both m values. We also show the fuU fields for both 
m values. Notably, for the case to = 1 the dipole mode is also excited, and eventually dominates the full field at late 
times. In the m = 2 case the smallest ^-mode that is excited is the octupole, such that the full field in the to = 2 case 
decays more rapidly than in the to = 1 case. We also show the local power indices in Fig. 1131 The late time values 
of the local power indices for each projected £ mode approach the same values independently of the value of to, in 
accordance with Eq. ([3]). 

In Table IV we re-organize the data appearing in Tables I, II, and III, so that the m-independence of the power-law 
indices is brought out in a sharper way. Specifically, the late-time decay rates for any allowed excited £' mode from 
an initial £ mode is independent of the value of to in all the cases we have examined. 
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FIG. 12: The full fields and their ^ = 3, 5 projections as functions of the time t for m — 1 (solid) and m — 2 (dashed). To help 
identify the modes, we denote the fields with e',m''pe.m- The initial data are for an initial pure octupole moment. 
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FIG. 13: The local power indices 3713 and 3725 as functions of the time t for m — 1,2. To help identify the modes, we denote 
the local power indices with i'^rnne.m- The case m = 1 is described by a solid curve, and the case m = 2 is described by the 
dashed curve. 



TABLE II: Scattering matrix for non-azimuthal modes with m — 1. Late-time power law indices (ine for the ^-projections of 
the fields for pure m = 1 mode initial data sets, £' = 2, 3, 4, 5. This table includes both "down" and "up" excitations. Asterisks 
relate to kinematically disallowed states. All figures are significant. 



Initial 


Projected 


Projected 


Projected 


Projected 


Projected 


£' mode 


e = 1 


£^2 


£^3 


£ = 4: 


£ = 5 


2 


+ 


6.997 


* 


8.995 


* 


3 


5.0001 




9.009 


* 


11.005 


4 


* 


7.0001 


* 


11.006 




5 


7.0006 




9.0001 


* 


13.008 
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TABLE III: Scattering matrix for non-azimuthal modes with m = 2. Late-time power law indices (ine for the ^-projections of 
the fields for pure m = 2 mode initial data sets, / = 2, 3, 4. This table includes both "down" and "up" excitations. Asterisks 
relate to kinematically disallowed states. All figures are significant. 



Initial 


Projected 


Projected 


Projected 


Projected 


£' mode 


e = 2 


£ = 3 


e = A 


1 = 5 


2 


7.0003 


* 


8.995 


* 


3 


* 


9.0006 


* 


11.001 


4 


7.0001 


* 


11.003 


* 



TABLE IV: Presentation of the data in Tables I, II, and III in a form that emphasizes the independence of the decay rate on m. 
Late-time power law indices e,mne,m for the ^-projections of the fields for pure m = 0, 1,2 mode initial data sets, £ = 2,3,4. 
This table includes both "down" and "up" excitations. Asterisks relate to kinematically disallowed states. All figures are 
significant. 



Initial 


Projected 


Projected 


Projected 


Projected 


£', m mode 


£ = 2 


£ = 3 


£ = 4 


£ = 5 


2,0 


7.009 


* 


9.01 


* 


2,1 


6.997 


* 


8.995 




2,2 


7.0003 


* 


8.995 


* 


3,0 


* 


9.009 


* 


11.008 


3,1 


* 


9.009 


* 


11.005 


3,2 




9.0006 


* 


11.001 


4,0 


7.002 


* 


11.008 


* 


4,1 


7.0001 


* 


11.006 


* 


4,2 


7.0001 


* 


11.003 


* 
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